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ABSTRACT 

A direct design method for solving the problem of 
robustness to cross-ccupling fperturbations in miltivariable 
control systems is presented. The method uses numerical 
optimization procedures to manipulate the system feedback 
gains as direct design variables. The manipulation is acccn- 
plished in a manner that produces desired performance by 
pole placement and robustness by modification of the mininun 
Singular values of tke system return difference matrix. 

Channels affected by cross-coupling perturpdation may be 
recognized by the character of their transfer function flot. 
The mechanism used Ey the pole placement and robustness 
routine in obtaining a robust design is evideat from the 
gain changes associated with the transfer function diagran 
and the zero shifts shown on pole-zero piots. The pole 
placement and robustness routine uSeS gain equalization and 
Zero assignment to modify the . characteristics of the systen 
in the areas of low Singular values, producing ea robust 
design. 

A modification of the jfole placement and robustness 
routine that may Ee applied to the design of rebust 
ebservers is also presented. Using feedback and filter gains 
as direct design variables a practical design proceaure for 


robustness recovery in observer based systems 1s obtained. 





TABLE OF CONTENTS 


i INTRODUCTION e e BS e s e a +d ] @ @ @ s = e s * s @ 12 
TT. Sligtie@NnrUl SINGLEWOUTPUT SYSTEMS . 2. «© © « «© « « 14 


rT . UAV ADE oY OTEMS e « 6 © 0« «© «2 © «0 © © #@ «© « 20 


Iv. ieee ees «6 6 « 4 6 6 © s 8 6 ss ew ttl lt ee SS 
Ve Cre eesAr TON, LESIGN FRCCEDURE “2. 3. « « « «© «© « «© « 32 
Vi. Pre eNOnt eer NOBLEM ls ss 4s «6 «© © «6 « s « « « » « 60 


Vil. Pe COrmaitnmmwoMmNb ELEY PROBLEM 5 < 2 s« «+ «ss « « « « 82 
BEBE OELS OBSERVER « « « s « « «© « © © «2 «© © © wo « 143 
IX. MmemUomeOesonVvEn DESEGN . s 2s « « « « « « «© « @« 124 
Xe CCNCLUSIONS .. ee 130 
Mmemeecr KErERENCES .« « « « © « «© ©» © © «© © @ ew ew ew ew) 618G 


mmm DLoOTREIBUTEON IIT5T . « « ss © «© © «© © ew ew wt lel w®€6 | 6UD 








———— 


Po SOF EASE so 


Cemparative Results Simple Problem. 


CH-46 Helicopter Parameter Definitions 


Design One Pole Placement 


Helicopter Prcblem Feedback Gains 


Okserver Parameter Data 





Biot OF FIGURES 


Basic Gain and Phase Margin Determination 

ot i iMiic «i's slic) «= ss <«  « .«. © 
Ge SoMiaclBGGGmP1Ot Gils = 96 s « « © « s -« 
lveitt=ctePlOteor stable System <«« « .« « « 
Mote very rerturpbed System 92. < .«. « « « « 
miGeavemNvaurst PlOE- Gs « .« « « « s « e 
Nyguist for Ineguality Additive Condition 

Minit oiredtave syYStem << = .« « s « « « « 
Nyguist Plot for Multiplicative System. . 
Typical Equivalent Feedback System... . 
PomdmpePioteot an Optimal System ..+-.. « 
Nga MCGDEOUE 9 se. es «6 < s « + s « 
Pacte miter —anput, Multi—output System. . 
NGiwnedive PCrEtuLbat2ON <0 5 5. « « « « « « 
MimEmmemca tive Per tirbDatwomus =< « s « < « « 
Design Space for Example Problem .... . 
Deceamme space £05 Column Probiem «. . « « « 
Mmuoctmarrye EXampie for Eteration . .. « 
Plot of Various Points With Zero Gradienot 

Sonstraimed Cptimization Example . . . « « 
Cmedieerzatron ol ADS Program . s. < « « « -« 
Comoenmsecction Dlagmram ss). s. « « « « s« + -« 
steepest Descent Algorithm . - . « « « « « 
Fletcher-Reeves Conjugate Direction 


DENG CIMACT es ses) 6 eRe « « «© «© *© « © « 


Universal Gain and Phase Singular Value Plot 


SMSCCEVeE I Melement ablOW  « « « «6 « « « -« 
Basten UlLtEerinput Ulta =-output System . . 


is, 
16 
a 
18 
19 
20 
21 
ae 
76 
24 
28 
XE, 
30 
g2 
34 
36 
38 
4 0 
G2 
Ly uy 
45 
48 


49 
27s. 
58 
61 





| 


a. 7 


Us 


URS 


Be 10 


ie 1 1 


Pola 
fe 13 


Mbeivanrraoye Nyauist £Or 2st3/(st1)2 .. .« 
Minimum Singular Value Plot, Example Problen 
ices ce OPAGkale Gh 1/St1 WM. . . 5 2 sw tl le 
Dome nicamovooeil 6 sli « =» « « © » « « « « e 
Singuiar Value Plot for Simple Problem. .. 
Clesed-loop Foles for Pole Placement Only 

Cc SCmmemere ells eis ls le e's) 6 es « «© + © « «© « 
pisiditeats Value Pllot Case Three . ..«+«+« « s 
Open-loop Transfer Function 2-1 for Baseline 
Transfer Function 2-1 Optimized( 3 Var) cue 
iMonsateCrerunetiOn I-11 Baseline . . .+«-e«-. - 
Transfer Function 1-1 Optimized (3 
WictNGUOMNC SiN s oe cellkclhetiel «<6 « «© « « «© e 
Closed~loop Fole-Zero Plot Case 3 .... . 
Pole-Zero Plot Simple Problem (4 Variables) 
Socom rock DlagEanm s/s 9s. << « « « « « «© e 
HeeanbaenmecONtrGleStrictUre 2 . < s 4 « « « e 
ES veccoumaceialneiMcMs ies) «6 » « © « « « « 
Alphatech Design Singuiar Value Plots ... 
Design One Perturbation Input ....+s+ « -« 
Destgueiwo Perturbation Input <<. . . « »« =» « 
Transfer Function $g- Jn Design One, 

NeMi@ mimes Ze Gs suciie 5 "sls = 0 oe « « « «© « 
Transfer Function 6a-d3 Design One , 

OMe nee CuMct ach cole MNSMESIEIG Tos « «4 «<0 « « 
Transfer Function éde- 58 Design One, 

Noi Oneninte ed eel Mslc lle clls o«0 « « « « ei <«.« 
Transfer Function $c- 39 Design One, 
Sceenieeo Um lec he aeus isi Mellel) « 6 6. cis 6 © «6 « -« 
POtcmeGOUr OES FOE SDESTOGH 1 . < « « « « « e 
Poca erOer Cres FOL Design 1 (Cont.) « « « « 
Hobe=GeborP#Ots fOr Design 1 (cont.) ... . 


Pole-Zero Plots for Design 1 (cont.) ... . 


93 


94 


Me) 


Sli 
98 
ss) 
100 
101 





1 | tials 
LP " 1 Lik . 








7.14 
7215 
ie 16 


ol? 
va18 
39 


72 0 
72 1 
8.1 
8.2 
8.3 
8.4 
8.5 
8.6 
a.7 
9.1 
9.2 
9.3 


The 68 to v Freguency Response, 


The S@to v Frequency Response, 


Transfer Function §¢-58 Design 
Nonopte. .« © « « «© « «© © © «© « 
Transfer Furction &- 5? Design 
Transfer Function &-ds Design 
Transfer Function &- da Design 
OMe zed «<2 6 « sss « « 
Singular Value Plot Design Two 
Time Response Design Two... 
Spaplhe Observer << .« « « « « « 


MprGMees tm INO 6. 2. lo eels) « «<6 


i 


1 


Nonopt. . . 
Optimized . 
Casez2, 


Case 2, Opt 


Two, Nonopt. 


Two, 


Nyguist Plot for Robustness Recovery ... 


Singular Values of Observer System... . .- 


Singular Valve Comparison Plot 


Nyguist for Computed Robustness Recovery .. 


Time Response Plot for Simple Observer .. . 


Observer Based Controller .. 


Singular Value Plot of Observer Results .. 


Time Response for Observer Systen 


102 
103 


105 
106 
108 


109 
110 
let 
es 
175 
117 
118 
120 
121 
122 
Pas) 
126 
ee 





ACKNOWLEDGEMENTS 


For over fifteen years the United States Navy has 
allowed me to make practical application of the science of 
flight as a Naval Flight Officer and as a flight test 
moo ject officer. I wish to express my appreciation to the 
Navy and the citizens of the United States for the cpportu- 
hity te further expand my knowledge of aviation by pursuing 
this degree. 

I would also like to express a Sincere tnank you to the 
Doctoral Committee of professors Platzer, Gawain, Latta, 
Strum and Collins for their support in these efforts. Of 
course special thanks must go to Professor D.J. Collins for 
heading the committee and having worked so long and hard 
with me cn two theses. Without his support this effort would 
never have reen completed. 

Thanks are also expressed to Professor G.N. Vanderplaats 
for allowing me to use the Automated Design Synthesis code. 
His help in integrating the code into the pole placement and 
robustness design program and the insight into the power of 
humerical oftimizaticn that he provided were extremely bene- 
ficial to the outcome of this thesis. 

The Naval Postgraduate School Library staf& and Roger 
Martin are appreciated tremendously for providing excellent 
service. No reguest for reference material was ever too 
Meeeicuit for them to handle. The W.R. Church Computer 
Center is acknowledged for their support inthe computa- 
tional processing for this thesis. I would also like to 
acknowledge the work of CDR Al Diel and LTs Cliff Cooksey, 
John Hoden and Mike Laptas for the work done to provide an 
excellent control system analysis pacxage at the Naval 
Postgraduate School. Their efforts made the analysis of data 
from the pcle placement and robustness routine much easier. 


10 





I would also like to thank my wife, Janice, and chil- 
dren, Bradlev and Melissa. AS anyone Knows who has under- 
taken an effort of this magnitude without their love and 
Support this would have been impossible. They have "gone it 
alone" guite often sce that this paper could be completed. 

Above all thanks to Him who started it ail. 


11 





I. INTRODUCTION 


With the rising interest 
brought on by 
to 


designer tc 


increasingly 


arisen develop design 
Specify systen 
time ensuring 


in 


relatively hi 


ness. the single-input 


PiehwiecemvdnlanlemcortEtol theory 


complex systems the need has 
methods that will allow the 
performance while at the same 


gh stability margins or robust- 


Single-output (SISO) case the 


designer has had the tools to do these tradeoffs in the forn 


of Nyguist, Bode and root locus plots. In the multi-input 


Multi-output (MIMO) case the classical methods are not 
totally appropriate. 

With the increased interest in MIMO systens numerous 
methcds of design have been employed to obtain suitable 


system performance and robustnesS with varying degrees of 


success. One primary method of design is to keep the ;flant 


as decoupled as possikfle throughout the design so that each 


individual element may be controlled independently and 
designed essentially asa single loop system. Rosenbrock 
{Ref. 1} has developed a procedure where the multilcep 
system is modified into a system that has diagonal 
elements that are much larger than adie OLs-cragond | 
€lements. This diagcnally dominant system is then in a forn 


where conventional Nyquist type techniyues can be employed 
MIMO design method is that 
method. This nethod 


optimization principles 


in the analysis. A third comaon 
of the 


euadratic 


uses a 
tO 


for various perfcrmance levels 


linear Quadratic (L¢) 
cost functional and 


allow the designer tc design 


cy adjusting the matrix weighting terms used in the cost 
Gane tion. The majer difficulty with all of the above 
methods is that they are not necessarily robust. This is 


especially true for cross-coupling terms between loops. 


AZ 








The primary achievement of this thesis has been the 
incorporaticn of tke time domain pole placement design 
procedure with a method of using return difference natrix 
Singular values to improve the robustness of the design. The 
technigue, which utilizes a medern optimization routine, can 
Significantly assist the designer in obtaining robustness in 
the face of cross-coufling perturbations. It has also been 
shown that the cross-coupling perturbation problem can be 
detected by using classical open-loop Bode diagrams as well 
aS modern control analysis. The pole placement and robust- 
ness design routine developed for this thesis has been used 
cn several problems discussed in recent literature. In these 
studies the pole placement and robustness design code has 
proven capable of meeting the desired goals cf pole 
placement and robustness and also brought to light scne 
interesting aspects of the cross-coupling perturbation 
problem. A slightly modified pole placement and robustness 
routine has proven effective in the design of robust 
cbservers. 

The remainder of tne thesis will present background 
material on SISO systems in Chapter Two and on MIMO systems 
in Chapter Three. Optimization will be discussed in Chapter 
Four aleng with a discussion of the Automated Design 
Synthesis (ADS) program used as the optimizer routine for 
the pele placement and robustness design procedure develored 
in this thesis. The thesis methcdology will be discussed and 
outlined in Chapter Five. This will be followed by charters 
discussing results from selected problems. Conclusions will 


re presented in the final chapter. 


23 





Il. SINGLE-INPYT SINGLE-OUTPUT SYSTEMS 

The purpose of this thesis research has been to develop 
a methoc of obtaining a robust multivariable control systen 
design. A brief review of the concept of robustness and 
stability in the framework of a conventional SISO systen 
will be dene before pursuing the concepts in a more 
complicated fashion in the following chapters. A simple 
interpretation of robustness is the ability of the system to 
tolerate design perturbations. These perturbations could be 
in the form of actuator failures, plant parameter uncer- 
tainty, unmodeled dynamics or nonlinear terms, or any one of 
Many other perturbations to the nominal design of the 
systen. : 

The primary reascn tor feedback systems is the control 
cf uncertainty within the systen. By appropriate use of 
feedrack, properties that would lead to an unstable systen 
may ke controlled. When stability and robustness asfects 
are considered for a SISO system, frequency domain design 
concepts, using either Nyquist or Bode plots, are normally 
used. Robustness in SISO systems is formulated naturally by 
the ccncept of gain and phase margins, both of which are 
readily availiable on the Nyquist or Bode diagran. 

In figure 2.1 a nominal feedback system can be seen with 
a perturtation element CO (s) placed in series with the 
nominal system. When of=1 the system is nominal and stable. 
To determine the positive phase margin of the system the 
value of @&(S)=QU(ju) = e?? vill be Changed dy varyingy 
until the system just becomes unstable. This value of Y will 
then Fe the system phase margin. The negative phase margin 
can be ccmputed in the same way. To find the gain margin the 


Magnitude value of & is increased until the system just 
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Figure 2.1 Basic Gain and Phase Margin Determination Model. 


becomes unstable. This is the upward gain margin. A sinilar 
dower Bargin is also defined. 

These gain and fhase margins may not be adequate meas- 
ures cf robustness [ Ref. 2] because they domnot aecoount for 
simultaneous variation in both gain and phase. Therefore, 
while large individual gain or phase changes may not desta- 
bilize the systen, small simultaneous changes in gain and 
phase may destabilize the system. This is not a major diffi-~ 
culty in classical SISO technigues because the effect can be 
easily detected. 

Gain and phase margin can be defined in terms of the 
cpen-lcor frequency domain plots in either the Bode or 
Nyquist format. Figure 2.2 depicts a classical Bode plot 
Showing gain and phase margin determination from the plot. 
The Nyquist plot may also be used to obtain this inforna- 
te OD . Nyguist criterion states that if th2 cpen-loop 
transfer function G(s)H(s) has n poles in the right half 
rlane and the limit cf G(s)H(s)=constant as s-»~2 then for a 


1s 





Phase~ deg Amplitude~ db 





Figure 2.2 Classical Bode Plot. 


stable system the locus of G(S)H(s) will encircle the -1+j0 
point n times in the counterclockwise direction as s varies 
along the Nyquist ccntour. If there are no poles in the 
right half s plane then the locus will not encircle the -1 + 
jO point. The diagram in figure 2.3 illustrates a neminally 
Stable system. The gain and phase margin may be determined 
directly from the diagran. 

Any change in the loop transfer function, provided the 
order of G(s)H(s) does not change, that changes the numter 
of times the locus of G(s)H (Ss) encircles the (-1,0) point in 
the Nyguist plot causes the system to become unstable. This 
leads to the conclusion that the minimum distance of the 
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PrEgure 92.3 Nyguist Plot of Stable Systen. 


locus of G(s)H(s) to the (-1,0) point 1S a measure of the 
system stability. ths distance concept carries over 
directly to the MIMC system as will be shown in the next 
chapter. Examples cf a multiplicative perturbaticn and an 
additive perturbation illustrate this idea. Figure 2.4 1s an 
additively rerturbed systen. Figure 2.5 snows the Nyquist 
meot for this systen. Assuming that the plant is itself 
Stable and the perturtations are also stable the diagram nay 
then be used to determine how near the system is to insta- 
bility for the given perturbation. 








Figure 2.4 Additively Perturbed Systea. 


Since the system is stable the (-1,0) point is encircied 
the correct number cf times by the nominal plant. If the 
locus of g{jw) in the diagram is warped until it passes 
keyond the (-1,0) point then clearly the number of encircie- 
ments of this point will change and the system will beccme 
unstable, assuming the order of the plant is not changed by 
the perturbation. To keep the locus of points from moving 
beyond the (-1,0) point eguation 2.1 must hold. 
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Figure 2.5 Additive Nyguist Plot. 
LAg(jw) Ii 1+ g (jw)! (2.1) 


mors) CONdition 1s illustrated in figure 2.6. The right-hand 
Side of equation 2.1 is just the magnitude of the return 
difference transfer function of the nominal systen. The 
multiplicative case is depicted in figure 2.7 with its asso- 
Cliated Nyquist plot in figure 2.8. The requirement for 
Stability is similar to the additive case and may be stated 
in eguation 2.2 


LAg(juw) 1 <i 1 + (9 (jw)7) 1] (2.2) 
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Pigure 2.6 Nyguist for Inequality Additive Condition. 


The aktove arguments will be applied again in Chapter 3 to 
develecp multivariable stability and robustness properties. 

The linear quadratic design has been the primary method 
employed in modern ccntrol design practice. In this method 
an optimal state feedkrack control law is develop2d to finda 
set of feedback gains that optimizes a chosen performance 
index. The performance index for the steady-state case is 
Seven in equation 2.3 
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Figure 2.7 Multiplicative Systen. 


5 all [= 


Po =-f xg x + u'R uj)dt (2525) 


where the x'x term and the u'u term form quadratics. The Q 
and R matrices are chosen by the designer to provide the 
rest comprogise between the minimum error of the system and 
the minimum energy needed to control the systen. The LQ 
method is ctased on the use of closed-loop state variable 
feedback fcr the control of the system. In the MIMO problen 
IQ methods have been used extensively because of their guar- 
anteed stability margins with diagonal weighting matrices. 
FOr diagonal weighting matrices the LQ method yields a 
guaranteéd phase margin of 60 degrees and -6 db to infinite 
gain targin. 
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Figure 2.8 Nyguist Plot for Multiplicative System. 
For SISO linear quadratic optimal regulators these 
Stability margins can be developed from inequality 2.4 


(Ref. 3] 


11 4 £'(jwI - a) mtkyj> 1 (2.4) 


We 





Writing the Kalman eguation as 2.5 


[eet £'(sI-a)-*bj|? = 1 + (42) 1 Gy(s) 1? (2.5) 
cne has further 
| 1+ G(s)H(s)1? = 1+ (1/,)16,(S) 1? (2.6) 


and for all s=jw and 0S w<eothe function (175) Gx (S) is 
greater than zero, therefore the Kalman inequality 1S shown 
to be 2.7. 


felt G(s)H(s)[ > 1 (Qe 





Figure 2.9 Typical Eguivalent Feedback Syste. 
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Mine graphical interpretation of this result is that the 
polar plot cf G(s)H(s) must remain outside the unit circle 


Semtered at the =1 + j0 point for all frequencies. Figure 


UNIT CIRCLE 





Figure 2.10 Folar Plot of an Optimal Systen. 


2-19 shows such a polar plot. Since the optimal regulator 
with the loop broken at the input to the plant (denoted by 
the x in figure 2.9 ) does not penetrate the unit disk about 


-1 + jO this means that the single input regulator will have 
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aphase margin of at least 60 degrees anda gain margin 
tolerance of fifty percent gain reduction and infinite 
upward margin. Further discussion of this property may be 
found in Anderson and Moore [Ref. 4]. Vinee dS Srasic 
review of the concerts of stability and robustness in the 
classical SISO system complete, the next chapter will extend 
some cf these basic ccncepts to the MIMO systen. 


Pays 











III. MULTIVARIABLE SYSTEMS 

Linear guadratic deSign has perforned relatively weli in 
aircraft control concepts because of the ability tc fernu- 
late the system State equation and to quantify acceptable 
perfcrmance indices for the system. Industrial applications 
for LC theory have been less successful and have led British 
researchers to look into forms of decoupled design methods. 
Cne of tke primary methods used in multivariable design is 
to make the system totally decoupled. This metnod allows 
each loop tc be designed as a separate entity by classical 
means. Cne of the primary difficulties with this nmethced is 
the proklem of finding a compensator which wili totally 
decouple the systen. The method also suffers ‘fron the 
effects cf£ cross-coupled perturbation terms. A method that 
does’ not tctally deccuple the system Dut makes the design 
robles simplier by designing a compensator that only causes 
the diagcnal terms of the transfer matrix to be dominant 
Omer all cff-diagonal terms has also been developed, 
[Ref. Dj. Classical frequency domain techniques are then 
used to design each loop of the system. The major difficuity 
with this and other single loop design technigues is their 
failure to account for cross-coupling perturbation terns 
that may interact between the loops. The individual loops 
Day be highly robust in these designs but the overall systen 
robustness may be lew because or the loop interacticn not 
accounted for in the design. This is the precise area that 
the singular value analysis has proven so beneficial in 
design . Singular value concepts may be applied tc conven- 
tionally designed systems to assess their robustness. for 
instance, if a system is designed by LO methods the designer 
May then formulate the transfer function of the system and 
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assess tke singular values of the return difference matrix 
of the system If the minimum Singular value is found to be 
low at somé critical freguencies the designer can then 
modify the Q and/or & matrices chosen in the LQ perforaance 
index and recalculate the design. In this iterative fasnion 
a rokust design would be developed. 

A generalization of the SISO Nyguist theory discussed in 
the previous chapter has been made for the MIMO proklen. 
This generalization leads directly to the singular value 
concept. The generalization is expressed in the form of the 
Multivariable Nyguist theorem which requires that a closed 
loop startle system have the same number of counterclockwise 
encirclements of the crigin by the locus of the det (I+tG (jeu) ) 
as the numter of open loop poles that are unstable. fThis 
theoren is formally stated as; 
let N{f(s)] denote tke number of clockwise encirclements of 
(-1,0) by the locus of f(s) as s traverses the contour D 
of figure 3.71 in a clockwise sense. The closed-loor systen 
will ke statle if and only if for all R sufficiently large 

N[f£(s)j)] = -P 
where | 
f(s) = -1 + det{ I+G(s)] = Yo(s)/ Ye.(s) -1 and 
P= the number of closed right-half plane zeres of 
Yo. (S) - 

The application cf the Nyguist theorem comes thrcugh the 
fact that a multivariable system will not b2 robust to 
modelling errors if the return difference matrix, eer is 
hearly singular for some frequency. Pa eG Ss eed. / 
Smmeuidr a small change in G may make I+ G exactly 
Seegudar.s This Causes the det({I + G) to become zero and the 
Nyguist encirclement count to change indicating an unstakle 
Meecha It iS possilie for very smail changes in I[+G to 
produce large changes in the determinant of I+G. The matrix 
I+G 


oa) 





Re(S) 





K: Zeros ell 


Paguiae ae0l Nyguist D Contour. 


10/sta 10/sta 
has dé€terminant 0.1/ (sta)?. If the element Pz, 15 changed by 


nae | 


only cne percent to 9.9/sta the determinant becomes 
1.1/7 (sta)? which is a significant change in the determinant 
wee. Ihereftore, it 1s evident that det (I + is net an 
accurate measure of how near the return difference is to 
Singularity. Researchers, in the Piers Of controls 
[Ref. 6], [Ref. 7], [Ref 8], [Ref. 9] have used Singular 
value analysis to determine how near the return difference 
matrix is te singularity. 

Since the number cf encirclements of the Nyquist diagran 
Changes as f(s) passes through the -1 point or when det (I+G) 
is zero it is imfortant to find how near the return differ- 
ence matrix I+G 1s to being singular. This nearness to 
singularity can be interpreted as closeness of the matrix G 
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Mmemene Critical point, -1. A quantity which can be used to 
express the nearness to singularity of the matrix 1s the 
minimum matrix singular value denoted by 9 . Given a matrix 


A the singular value may be expressed Dy equation 3.1 
. a. J 
C (a) = min ( ); (A A)) (ae) 
- 


where mw (ATA) is the eigenvalue of the complex ccnjugate 


transpose of A times 3. A basic MIMO linear system is 





Figure 3.2 Basic Multi-input, Multi-output Systen. 


Gepacted in figure 3.2. An additive perturbation to the 
plant is shown in figure 3.3. If the plant is stakle before 
the perturbation is added to the system the Nyguist theoren 
will te satisfied and the locus of GH will not encircle the 
omeo Critical*point. When the perturbation is added to the 
System as long as the Nyguist locus is not forced to 


Face, 





t) 





Figure 3.3 Additive Perturbation. 


encircle the -1 point the system will remain stable. A 
sufficient condition, recalling the SISO discussion in 
chapter 2, for the perturbed Nyguist plot not to change 
encirclements is that the norm of the perturbation AG 


remain less than the norm of the return difference matrix as 
expressed in equation 3.2. 


ieee G (Ge) 11< 1711 (1*G)~*) 11 (3.2) 

ww 20 This conditicn will guarantee that the Locus of the 
det (I+G) does not pass through the -1 point. If the l, or 
Euclidean norm is assumed for this condition the equation 
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3.2 may Le expressed in terms of singular values as equation 
ea « 


o(AG)S C(I + G) (3.3) 


This result states that as long aS the nmaximim singular 
value of the perturkation matrix AG is below the nmininua 
norm value cf the return difference matrix the system will 
remain stable. The problem of guaranteeing robustness 
becomes that of finding the largest norm of the perturbation 
guantity, the largest singular value, for which the smallest 
norm or singular value of the return difference matrix will 
remain ncn Singular. 

The multiplicative form for a system such as figure 3.4 


gives the similar norm eguation in 3.4 


Li O@(je) bis vb (ZL + (G)- 4) 721 (3.4) 
CJ 20. which may te expressed as; 


Set Ac) < © (1 + G-! ) (25) 
Singular value decomposition software is readily availatle 


to détermine how near the tmatrix I+G or I+ (G)-? pie xe) 
Singularity. 
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| Figure 3.4 Multiplicative Perturbation. 
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IV. OPTIMIZATION 
The yurrose of this chapter will be to briefly describe 
several of the currently employed optimization technigues 
and the Autcmated Design Synthesis (ADS) program, [Ref. 10], 
which emyloys these techniques. In general, optimization 
implies finding the "best" possible solution to a frotklen. 
In actuality the best solution found by an optimization 
technigue might really only be a "better" solution to the 
problem. The purpose of ADS and other optimization routines 
as to allow a rational search to be conducted to find the 
Fest possible design. Tne techniques of numerical optiniza- 
tion are used to logically vary the various parameters that 
affect the design until a good solution is found. 

As an example of an unconstrained optimization prcblen 
consider the following problem developed in [ Ref. 11]. The 
problem is to minimize the function 


eee 10x = 20K, 2X, +10K 2 +X ,2- 2x, +5 (4.1) 


fry is ctten called the objective function, the cost func- 
tion or the penalty function. Since there are no conditicns 
imposed on the design variables, x, and x,, and no addi- 
tional ligits imposed on the overall design, the pfroblem is 
considered to be one of unconstrained minimization. Figure 
4.1 represents this problem in the design space. From the 
figure it appears the optimum is near the point i,1. 
Calculus may be aprlied to determine the optimum exactly. 
Taking the derivatives 
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Pigure 4.1 Design Space for Example Problen. 


OF (2) /)x, =40x, 3-40 x, x, +2x, -2=0 (4.2) 


OF (x) /Ox,=-20x,2 +20x,=0 (4.3) 


and then sclving the set of equations it is found that x, 
=1.0 and x. =1.0. 

If design conditicas are imposed on the problem then the 
optimization becomes one of constrained function nmininiza- 
tion. In other words, while the minimum of a function is 
Berti sought, this finimum must exist within the limits 
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imposed ky the design conditions or constraints. The tinitua 
of an unconstrained function will not necessarily be the 
sane Dinimum for a censtrained function. 

As another example , also from [Ref. 11], a design is 
sought which gives the minimum weight of a particuiar 


column. This weight is expressed as; 
W= pah= e/tDh (4.4) 


where ie is the unit weight of the material and A the cross 


sectional area. The stress in the column is given by: 

O =P/A = PSR Dt (4.5) 
without coing into detail the design will be constrained by 
the allowarle stress cn the structure. Other constraints , 
Euler buckling and shell bucxling, are also of interest to 
the designer of this column. The design probiem is then 
stated as: | 


Minimize W = 7 Dth (io) 


mene the constraints 


g(1) =" yo- 1s 0 (7) 
g({z) =S 4-15 0 (4.8) 
pe) = 97G, ~ 15°70 (4.9) 


With D2 10e-06 and t 210e-06. 
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g(4) = t-Ds 0 (4.10) 


om 104°C, «wefer to the allowable design stress, fFuler 
buckling stress and the shell buckling stress respectively. 
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Wall thickness rario, tho 


Figure 4.2 Design Space for Column Problen. 


Figure 4.2 is a design space diagram for this problen. It 
is noted from the figure 4.2 that the optimum is not unigue 
and can ke any value along the arc noted as optimum in the 
figure. 
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The general optizization problem is then written as: 


Minimize F(x) objective (Qa 7) 


sukject to g (x) s OQ j=1,m inequality constraint (4. 12) 


hy (x) = 0 k=i1,1 equality constraint (4213) 


x. S SE Ac ay i=1,n side constraint (4. 74) 


Seeeeea = COL (X1,X2,-..,%N) 1S the design variable. 

The methods used to solve this problem are usually 
Gterative. After the establishment of an initial set of 
Variatles the optimizer will update this initial value until 
the cptimum values are found. Again ,borrowing from 
{Ref. 11], the iterative technique may be demonstrated by a 
Simple exarple. Figure 4.3 is used to illustrate this 
peoclem. Given the initial cata set £° the formula 


me) yO se o* s1 (4. 15) 


can be used to upgrade this estimate of X. The vector S is 
the search direction for the iteration and the scalar quan- 
tity ® is the distance of the move in the S directicn. 
fee g at X° at a= Gesired to reduce the objective func- 
tion. The search for values of X that reduce the objective 
Bimetion 1S made in the § direction which in this example is 
PBiewOrrosite of the gradient of the function at point X9. 
im choice of S could be arbitrary as long as it reduces the 
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Figure 4.3 Illustrative Example for Iteration. 


value of the objective. If s is chosen to be opposite the 
gradient of the objective then the search would te a 


steepest descent search. Once the direction is chosen then 
the scalar Must be found to minimize the objective alceng 
this direction vector without crossing a constraint 


Poundary. The technique for finding & is to evaluate X, the 
objective and the constraint functions for values of & using 
computer code and numerical interpolation to estimate x. 
This & value then gives the pinimum value for F(X) along 
this s¢€arch vector. Now that «x for this direction is known 
the new eyuation beccnes: 


x2 = x1 + ose (4. 16) 


ema mow X! is used to start a new search in directicn §* 


and compute a new a*that continues to reduce the objective 
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function. These one dimensicnal searches are continued until 
no more reduction in the objective can be found. hot o2s 
point no further design improvement is possible. 

Cetinization techniques do not always lead to the abso- 
lute cptimuz when applied to problems of practical interest. 
The reasons for this could be numerical ill-conditioning of 
the problem formulation or simply that there are multifle 
solutions tc the problem. Because of these difficulties it 
may be adviseable to choose several different starting 
points for the optinmizaticn process and use engineering 
judgement as to the design most applicable to the preblen 
under analvsis. 

Considering the unconstrained case first where the 
desire is tc minimize the function F(x), it is well known 
Pec F(x) Willi have a Binimum where the gradient of F(x) is 


Zero. That is: 

grad (F (x) ) =0 (4.17) 
with the Y F(x) defined as; 

mead {F(x))= ( OP (x)/oxX,,-+-, OF(X) /ox,) ' (4. 18) 


Figure 4.4 shows why this is a necessary condition but dces 
not guarantee a glokal mininmun. the gradient of F(x) is 
zero at all four points A,B,C,and D. dowever, only A and D 
are Minima. A would re the global minimum for the function 
as defined here. D would Fe only a relative minimum. To 
check that the zero gradient corresponds to a minimur the 
Hessian matrix,i.e the matrix of second partial derivatives, 
can be examined for fesitive definiteness. A positive defi- 
hite Hessian ensures a relative minimum. The only way to 
rrove a glokal minimum for the function is to show that the 
Hessian matrix is positive definite for all design variables 


X. A test that is seldcm possible to perforn. 
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Figure 4.4 Plot cf Various Points With Zero Gradient. 


If the problem of minimization is a constrained problen 
the situaticn is different from that discussed abcve. fthe 
objective function gradient does not have to be zero at the 
optimum. Figure 4.5 illustrates this case. Using figure 4.5 
and assuting a start at point A it is necessary to choose a 
Search direction S which will reduce the objective function 
while not violating the constraint functions. Any direction 
that will reduce the objective is said to be useable. This 
is seen to be the half plane sector below the tangent to 
F(x) at pcint A. If a half plane to the right of the 
tangent to the active constraint at point A is ccnsidered 


then the feasible sector is formed. A combination of the two 
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conditions gives the useable, feasible sector where the 
search direction S must be chosen. In mathematical nomencla- 


ture the above argument may be stated as; 


uscable direction grad({(F(X))eSs 0 (4. 19) 


feasirle direction grad ( g,(x))° SE WM petepe eulah 5 (4.20) 


for which J; (x) = 0. 

Point B in figure 4.5 shows a point where the gradient 
Sete OL JEective and constraint point in exactly the oppo- 
Site direction. At a point such as B the only search vector 
S that meets reguirements for useability and feasibility is 
tangent to the constraint boundary and to a line of constant 


ebjective function. This condition is stated as: 


grad(F(x))+2@Aygrad(g;(x))+2r grad(h,(x)) = 0 (4. 21) 


where A4;> 0 and Ags, unrestricted. 

With this brief kackground in optimization complete the 
Automated Design Synthesis (ADS) program will be briefly 
discussed. More detail on the ADS routine can be found in 
[Ref. 10}. This code was developed as a follow-cn tc the 
successful CONMIN code [Ref. 12] developed by Vanderplaats. 
It 1S designed as a black box optimizer which allows the 
user tec chcose combinations of one dimensional search, opti- 
mization algorithm and optimization strategy. These will be 
discussed later. For the user with specific requirements the 
code may be tailored by parameter modification to meet 
Specific requirements. for the work done in this thesis it 
was assumed the user of the code has no detailed knowledge 


of oftirization and will want to use the code in the 
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Figure 4.5 Constrained Optimization Exaaple. 


Simplest form. As such, the ADS calls made from the main 
program use only defatlt parameters and first focward finite 
difference gradients. Should analytical gradients be avail- 
able the user could use them within the code if desired with 
moe difficulty. The calls made by the user to the ADS 
routines specify several important aspects of the preblen 
solution. The user may call for an optimization strategy to 
ke used in the routine. This is not required and its use 
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depends cn the problem. These strategies are discussed in 
[Ref. 10]. Two strategies used most often in the analysis 
done for this thesis have been sequential unconstrained 
Minifization using quadratic exterior penalty function and 
the augmented Lagrange multiplier (ALM) method. Others 
available include seguential linear programming and seguen- 
tial guadratic programming. 

The kasic optimizer is also chosen by the user from two 
unconstrained and three constrained optimization algcrithms. 
The unccnstrained algorithms are Fletcher-Reeves conjugate 
directions, Davidon-Fletcher-Powell (DFP) variable metric 
method and the Broydcen-Fletcher-Goldfarb-Shanno (BFGS) vari- 
able metric methced. The method of feasible directions and 
robust feasible directions are available for constrained 
Minimization. 

The user has available several types of one dimensional 
searches using Golden Section orf polynomial approximation 
technigues. The ADS code has tailored these one dimensional 
search algcrithas fcr the unconstrained and constrained 
cases, allowing the user to make appropriate choices fer the 
type of problem to re solved. Figure 4.6 shows the kasic 
organization of the ALS progran. 

Since it is not the purpose of this chapter to ccver the 
subject of optimization in detail all the possible routines 
in ADS will not be discussed. A few of the routines found to 
work well for the werk completed in this thesis will be 
briefly described, however. There are several metacds of 
optimizing functions cf one variable or one dimensional 
searches. For instance a large number of points, n, could be 
chosen and the functicn F(x) evaluated at each point. The 
point corresponding to the smallest value of F(x) could then 
ke considered the oftimum value of the one dimensiozal 


search. This method is hit and miss and better methods ope 
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Figure 4.6 Organization of ADS Prograaz. 


locating optimum points can be used. The methods used in ADS 
are Golden Section ard polynomial interpolation. 

Golden Section search methods are easy to program on the 
digital computer and do not require continuous derivatives. 
They kave a known ccnvergence rate and are reliable for 
poorly conditioned froblems. The major drawback of the 
Golden Section routines is the large number of function 
€valuaticns require. The Golden Section algorithm is simfly 
illustrated through the use of figure 4.7 Assume that X9 and 
Xi are known to be kounds on the curve's minimum value. 
Also, the function values of F(x), F° and F! ,are evaluated 
and known at these foints. By picking two intermediate 
points X2 and X3 where X2 < X3 and evaluating the function 
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Pigure 4.7 Golden Section Diagran. 


at these pcints a fccund on the minimum may be modified. 
Since F2 for this figure is larger than F3 the point X? 
forms the new lower bcund. The minimum is now between X2 and 
Xi . If the function F* at X* 1s determined and shown to be 
larger than F2% then X* becomes the new upper bound on the 
Dininuna. By repeating this procedure the bounds may be 
Narrowed to any desired tolerance. The Golden Section rule 
is applied to this problem to reduce the bounds in the 
quickest possible time. By appropriately picking the values 
of X's at which each function evaluation is made an effi- 
client algorithm that uses the ratio proportion of the Gclden 
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secticn,i.e. X*/X! = 1.62803 is developed. By choosing a 
meee 2 based on this Golden Section rule where @ = 0.38197 


estimates for interior points X2@ and X3 can be made as: 


me= (1 - £0) xo +@x1 (4. 22) 


XxX 3 


aero + (1 = @j)x! (4.23) 


As each new bound is found the process repeats until the 
accepted level of convergence is reached. 

The polynomial interpolation method 1S accomplished by 
first fittirg a polynomial curve to the points about where 
the mininum is desired and then finding the minimum of the 
polyncmial function. For example if the function F(X) is 


approximated by a guadratic as: 
X2 (4.24) 


Then the value of xk, X , where Ft ais zero can te shcewn 
Hobe : 


x* = wa, /2a, (4.25) 


He a 


polyncmials may be used in Similar fashion. 


5 is positive then F will be minimum. Other degrees of 


Now that the basic one dimensional search methods have 
been reviewed the next step is to examine the basic cpti- 
mization routines. First, the unconstrained case will be 


reviewed. The optinum x*is at the point where 
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grad( F(X ))= 0 (4. 26) 


Several zere order methods exist for the purpose of finding 
the mininum value. These include random search, Poweli's 
method and’ Box's' method. Since ADS does not use these 
methods they will not be discussed. 

Autcmated Design Synthesis makes use of first crder 
methods which will new be discussed. The steepest descent 
method is best known, but poor in performance. Just as the 
Mame implies a search direction is chosen opposite the 
gradient of the objective function. At iteration 0 


s° = = grad{ F(X9%)) (4.27) 


Figure 4.8 shows this algorithm geometrically. Note that 
the methcd simply stair steps its way down the "hill" to the 
valley or minimum. ALS uses the Fletcher-Reeves modification 
to steepest descent. In this routine a conjugate direction 
is chcsén to improve the speed of the search. Figure 4.9 
Shows how this method tracks to the solution. 

The variable metric methods listed earlier are usually 
more frowerful than Fléetcher-Reeves because they stcre infer- 
mation that allows tke algorithm to approximate the inverse 
of the Hessian matrix or second derivative. For further 
discussicn cf these methods see [Ref. 11] or other similar 
optimization texts. 

ATS employs two direct methods for constrained minimiza- 
tion. One method is that of feasible direction and the cther 
is the methcd of robtst feasible direction. Since these two 
methcds were seldom employed within the work presented in 
this paper they will not be discussed. 

The methods chosen to handle the constrained minimiza- 
tion proklems formulated in this thesis are referred to as 
Sequential Unconstrained Minimization Techniques (SUMT). 
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Figure 4.8 Stee pest Descent Algoritha. 


SUMT methods are methods which formulate the okjective 
function and the constraint functions into an augmented 
objective function and then solving the problen as if it 
were an unconstrained optimization tasks ADS employs 
interior and exterior penalty function techniques as well as 
an Augmented Lagrange Multiplier (ALM) technigue. 

The exterior penalty function method is incorporated by 
fOrMing a penalty from the constraint equations. This 
Penalty is cf the forn ; 


{ 
P(x) =F max (0,g, (x))2#=& (hy (x)) 2 (4. 28) 
Jat Ke} 
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Figure 4.9 Fletcher-Reeves Conjugate Direction Algoritha. 


From this equation it is obvious that P(x) is zero if all 
the ccnstraints are satisfied. That is, if g; (x) §$ QO and h, 
(x) = 0, then all conditions of the penalty function are 
Satisfied. If an element of the penalty function is violated 
then the penalty increases as the square of the violated 
constraint. A pseude or augmented objective function is 


formulated where : 
PlXety)= FH) + Cp P(x) (4. 29) 
The ccnstant ©, is a weighting factor for the penalty. It is 


adjusted with ADS as the optimization proceeds to allow the 
program to systematically converge to an optimum solution. 


49 





One disadvantage of this method is that if the optimizer is 
stopped short of the optimum the design will be in the 
infeasible region and rrokably not useful. 

Criginally developed as a method to solive eéguality 
constrained problems the version of ALM in ADS has Deen 
modified to work with both equality and inejguality 


constraints. The statement is as follows: 


minimize F(x) (4. 30) 


sukject to hy (x) = 0 k=1,1 (4.31) 
Next a Lagrangian is created such that 
g ~ 
L(x) = F(x) + ZA. LBS (4. 32) 


Then the eguation is augmented with an exterior penalty 


function such that 
¥¢ 
A(X, A 05) F (X) * 214,b, (OSE tte 5 Gane) <a) (47.33) 


The power of this method is that in thedry precise 
constraint matching is possible whereas inthe exterior 
penalty function metnced it is not. The full details of the 
methcd will not be ccvered here, however, the final ferm of 


the ok-ective functicn will be included: 
yr Zz i 2 
A(x,A,Tp )=F (2) +2 Osimy)+ 2 bela) trp (hy (ey ) (4. 34) 


where ¥.=max(g 5 (x), -5/20p ), and the update formulas for 
the Lagrange multipliers are: 
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Pi? 
d; =) i 2c, (max (g ; (xX) ,-Aj/2r,) ) Je pe (535) 


Pri r - 
ae =) + 2rphy(x)  k=1,1 (4. 36) 


With this Erief review of optimization and the ADS progran 


completed the next chapter will discuss the program develop- 
ment for this thesis research. 
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V. OPTIMIZATION DESIGN PROCEDURE 

Tc accomplish tke task of designing a control systen 
with acceptable time domain performance and rotustness 
characteristics in a straight forward manner i numerical 
optimization technique was chosen as the method of ismplemen- 
tation of the design algorithm. Using numerical procedures 
to adjust selected design variables, in this case the feed- 
tack andyor filter gains, a desired level of performance can 
ke achieved. This level of performance is actually a ccembi- 
nation cf time dcmain performance and robustness or 
frequency dcmain performance. By establishing the criteria 
for the system perfcrmance in terms of an optimization 
objective and constraint functions a versatile procedure can 
ke developed to set the system feedback gains and affect an 
acceptable design in terms of performance and robustness. 
The pc-.e placement and robustness (POPLAR) design procedure 
uses pole placement to establish a designer selected 
perfcrmance level and then a minimum singular value level to 
estaklish the robustness. 

The pole placement portion of the procedure will be 
discussed first. The pole placement technigue was chosen 
because it was relatively easy to implement through a numer- 
meat Optimization routine. By using this numerical procedure 
it is also simple to incorporate robustness into the proce- 
dure along with the ferformance reguirements. A numerical 
technigue similar to one posed in [Ref. 13] was chcesen for 
the pele placement algorithm. An unconstrained optimization 
routine frem the IBM IMSL library was used for this rprogran. 
The rcutine, a Newton method, uses adjustments to the output 
feedback gains to reduce the size of an objective functicn. 
This cbhjective functicn was expressed as a function of the 
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pole locaticn of the system and as the objective was reduced 
the yoles were moved toward the desired locations. TO 
provide more versatility in the pole placement algorithm a 
methed that can use constraint functions as well as uncon- 
strained optimizaticn was chosen for this procran. The 
designer may use eitker objective, constraint or a comrbina- 
tion cf functions to secure the desired pole locations. As 
currently implemented in the program the cost or objective 
porticn of the pole placement procedure is constructed as 
equation 5.1 


A 
OBJ = Z (dg. - dg; 2 + (Ary, -Az; )? Boy) 
e=/ a 
where Ax = real €igenvalue 
»; = imaginary eigenvalue 
M, = desired eigenvalue location 


Mr, = desired eigenvalue location 
Be Constraint formulation is a function taat must be 
kept negative or the constraint is violated. It is written 
as equation 5.2 


(3) = V Oe 7 dey)? + Card)? 7 © (5.2) 


where r is a tolerance circle established as a function of 
Pele placement position. Since the aim of the optimizer is 
to keep g negative any time the » function of the 
constraint is greater than r the constraint will beccme 
active, i.¢€. violated. The optimizer will then attempt to 
move the constraint tc the inactive status by adjusting the 
design parameters of the systen. 

Consideration of implementation of the freguency demain 
or rokustness porticrn of the design procedure begins with 
the concept of MIMO phase and gain margins. Several useful 
theorems on singular value analysis of multiloop systems are 





presented in[Ref. 8]. One of these theorems relates the 
Matrix Singular value of the return difference functicn tc a 
parameter,®, and further shows that as long as the taximun 
Singular value of the perturbation function (L7~! - f) 
remains less than this «, the system remains stable. The 
value of & is then related to gain and phase margins of the 
MIMO system. The relationship developed is given in egua- 
tions 5.3 and 5.4: 


gain margin = GM =1/ (1+&o) (Sx } 


phase margin = PM = £tcos 7! (1-&,2/2) (5:64) 
provided that equaticn 5.5 holds. 


G(i+G)2 Ko (5 25) 
fer scme “X,<1 
These phase and gain marginS are guaranteed in every loop 
Simultaneously. 
Universal gain and phase margin curves, [Ref.j 14], tased 
on the minimum singular values of the return difference 


matrix are developed from eguation 5.6. 


(1-1-1) =nax (1- 1/k,) 2+ 2/k, (1-cos ¥n) (5.6) 


Mer all n with kK, > OQ. These curves shown in figure 5.1 
allow the deSigner tc pick a singuiar value that corresrfonds 
to a specific gain and phase margin for a given systen. In 
addition tothe universal gain and phase flot [Ref. 15] 
develores an optimizer solution for formulating a robust 
contrcller using the CONMIN algorithm [Ref. 12]. 
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Figure 5.1 Universal Gain and Phase Singular Value Plot. 


Since the universal curve in figure 5.1 provides a 
convenient method of specifying gain and phase margins in 
terms of Singular values the robustness portion of the pole 
placement and robustness design procedure uses the mininun 
Singular value level of the return difference matrix to 
determine the robustness. The minimum Singular value level 
is formulated as an objective or constraint function in 
equation 5.7 


J= 2. (max (0, (S,- Co (jo,P)) ) ? (5.7) 
=) 
The optimization procedure may be used to change feedback 
gains until the Binimum Singular value is raised above this 
desired design level. Although the same formulation can be 
used aS a hegative constraint function it has not’ been 


implemented as such within this program. There are numerous 
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ways the singular value formulation could be iznplemented 
within the rfrogram by changes of the code if design xequire- 
ments forced such changes. 

The pcle placement and robustness design program is 
based on the ADS code discussed in Chapter 4 to isplement 
the désign variable selection procedures. The pole placement 
and robustness program consists of two separate programs. 
The first program is used to provide designs for state or 
output feedback problems while the second program is’ used 
for okserver or filter designs. 

For the state or cutput feedback design program the user 
must input the plant matrices A,B,¢€ and initial starting 
values for the feedback matrix F. The matrices correspond to 
the fcllcwing linear differential systen: 


x= Ax + Bu (5-8) 
ei. (S300) 


A feed-fcrward matrix has not been considered in the current 
program. A feed forward matrix could be added td the froce- 
dure if required for specific design cases. 

As the design prcgram is currently coded the user may 
run ouvutput feedback or state feedback design by specifying 
the C matrix as the diagonal(I) matrix for state feedtack. 
The program relies on initial starting values of the feed- 
Fack gains, fF. As discussed in Chapter 4 there is no 
guarantee that the oftimum found by the procedure each time 


is the glokal optimum or that the procedure will always 





converge to an accertable solution. The ability tc select 
acceptable starting values for the feedback gains will make 
the procedure more efficient in operation. MomCuL cit) / 
employed, the prograz is used to obtain pole placement and 
robustness for a given set of Starting gains and a selected 
optimization routine from the ADS program. If the optimizer 
is net atle to meet the desired design goals on this vrogranm 
run two cptions are available. First, change to a different 
optimization routine from the list of available ADS routines 
and rerun the problen. This was usually successful in 
imprcving the design. Second, the designer uses 3 new set of 
Starting values for the feedback gains and repeats’ the 
design procedure. Both options might be used on particularly 
Stetacult cases. 

The pole piacement and robustness design procedure has 
consistently been able to find improved designs; however the 
program does not always yield acceptable designs. Certain 
problems require changes in the optimizer routine and modi- 
fication in the initial feedback gain starting values in 
order to obtain acceptable designs. Using the IBM 3033 time 
Share system the ele placement and robustness routine 
requires akout ten CPU seconds to work a second order 
BeoolLem and on the crder of 15 to 60 CPU seconds to run a 
forth order problem. The actual amount of time varies with 
optimizaticn requirements and time share utilization. 

The observer rcbkustness design program requires two 
passes of the ADS pregram. In the first pass the feedback 
Sains, fF, of the controller are computed to ottain the 
desired role locations. The second pass of the ADS routine 
1s used to adjust the observer gains to recover the systen 
robustness. The twe pass procedure was chosen because it 
allows a smaller numter of deSign variables at each stage of 
the ortimization and much more efficient computer usage. 
Figure 5.2 shows how the observer is implenented. This 


diagram is algebraicly stated as; 


57 





x =fA “a xl + i aScetol 
g - Ae-BK A B 
A. = A - FC 
maerbesx 1S the state, x IS the estitiator variabl>, F is the 
feedback gain and K is the observer gain. The design cf the 


feedback gains and the observer gains are accomplished as 


Separate guantities in keeping with the separation frin- 





X=AX+KY + BU 





Figure 5.2 Observer Implementation. 


Ciple. In using the pole placement and robustness design 
procedure for the observer systen, initial values of the F 
and K matrices must re input. The same or different opti- 
mization technigues from ADS may be enployed. 

The routines contained in this thesis have been Eased on 
the input additive singular value level. The pole flacement 
and robustness deSign algorithm computes input additive, 
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output additive, input multiplicative, and output multipli- 
cative singular values. Any of these singular values can be 
incorrorated into tke objective or constraint formulaticns 
but have net been for this version of the progran. 

In summary, the pole placement and robustness design 
procedure is a straight forward numerical optimization 
procedure for the practical application nodern MIMO systen 
analysis. The new asfects of the procedure are the implemen- 
tation of both pole placement and robustness criteria within 
the same design program. The versatility of the pole flace- 
ment and robustness design is obtained Dy incorporating a 
state of the art cftimizer routine ADS, with currently 
available singular value computation routines. Using the 
cptimizer format fcr the pole placement and robustness 
design gives the designer the ability to modify variables 
directly that affect toth the time domain or performance of 
the system and the frequency domain oor robustness of the 
system. The numerical optimization incorporated into the 
pole placement and srcbustness design is flexible enough to 
incorforate the double pass design technigue using the sefpa- 
raticn of the feedback and filter gains to obtain robustness 
recovery for cbserver based controllers. 
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The purpose of this introductory problem is to review 
some results of modern multivariable robustness theory using 
a Simple problem. This same simple system is then used asa 
test problem for the optimization technique develored for 
this thesis. The protiem provides excellent insight into the 
cross-coupling problem and demonstrates how effectively the 
pole placement and i«Icbustness design procedure can be in 
dealing with the cross perturtation terms. 

The proklem chosen for this introductory analysis comes 
mpom f[ Ref. 9]. Paguner O.1 Sed diagram OL this Easic 
systen. In this proklem a simple plant is specified by the 
following linear systen: 


ep |~1 Ope +11 b,,] fo (6.1) 
x2 0 -1 |ix2 0 u2 

where : 
ai = xt (6.2) 
y2 = X2 


A ‘fe¢eedkack compensation of the form of eguation 6.3 was 
assumed. 


mi = =~ | Xt} + juct (63) 
i 2a A We 


Which gives a closed-loop system, equation 6.4 


x1] = [-2 Bee aml [ue (6.4) 
x2 0-2 ree 0 1 uc2 


This system has eicéenvalues at -2 ,-2 andis therefore 
Stable. Using equaticn 6.5 
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G= C(sI-2)-*B (6.5) 
the transfer matrix may be written as equation 6.6. 


G =]1/st1 b,-/s+1 (Gao) 
0 1/st+1 
which gives the return difference matrix, equation 6./7. 


I+G(s) =|st2/sti [,, /st1 (6.7) 
0 ¢+2/s+1 


Y2 





Figure 6.1 Basic Multi-input Multi-output Systen. 


The problem shows the inadequacies of classical methods 
in establishing the robustness of the system. A brie? review 
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cf the results will kre presented. Using the return differ- 
ence matrix (1I+G), the determinant may be written as 


equation 6.8. 


eg 1G) -1=2s5+3/ (s+1) = (6.8) 


A 


Therefore, the multivariable Nyguist diagram will be as 


Im 
_. _ - ee 
i mG 3 Re 


Figure 6.2 Multivariable Nyguist for 2s+3/(st1) 2. 


shown in figure 6.2. The diagram does not encircle the 
f= 1,0) point and is indicative of a closed-loop stable 
system. Considered as a SISO system one has a gain margin of 
-1/3 to @ and a phase margin of +106 degrees. Under this 
Criteria one can conclude that this 1S a good design. This 
will te shewn later nct to be the case. 

Multi-input systems are often designed a loop at a tine. 
Applying loop design to the system the transfer functicn for 
either lcop taken sefarately becomes equation 6.9 
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G=1/s+1 (6.9) 


Figure 6.4 shows the Nyquist diagram of this problem. Using 
this Nyguist diagram the system is indicated to be stakle 
and have fhase and gain Margins of GM=(-1,@) and 
FM=(+180). This analysis does not show the true nature of 
the robustness of the system. Since the factor b 1S Tot a 
parameter in either of the Nyguist curves it plays no part 
in the stability determination using these diagrams. 

Using the criteria of singular values discussed in 
Chapter 3 a measure cf the nearness to instability for this 
problem may be obtained by plotting the minimum Singular 
value of (I+G),@. (For numerical calculations a value of 
B,g = 50 is assumed). Figure 6.3 shows the plot of this 
value vs. freguency. This gives a minimum singular value of 
about -23 db or 0.071 which corresponds to a gain margin of 
about 0.93 to 1.08 and a phase margin of +4.1 degrees. These 
phase and gain Margins are guite small and are evident in 
the cress-feed perturkation problen. 

A perturbed system as Shown in figure 5.5 can be 
Froduced which will lead to stability problems with very 
small values of perturbation. The closed-loop systen if 
perturked by a smali perturbation, o/b, , where by, isa 
large number, will have as a characteristic equation, eéqua- 
tion €. 10 


(sI-A)= s2 +¢ Us +6 (6. 10) 


with the eigenvalues of s=-2+ ‘1/5. There is one positive 
Toot in this solution and the system is unstable. 

To determine the nature of the robustness of this system 
the return difference matrix must be considered. If the 


return difference matrix of the transfer function (I+G (ju)) 
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Figure 6.3 Minigaugz Singular Value Plot, Example Problen. 


1s nearly singular at some frequency @, then the multivari- 
able syster will not be robust with regard to any modelling 
errors within the systen. This is because any small change 
in G(jw) can then make I+G(jw) singular and the det (I+G) 
becomes zero, thus changing the encirclements of the Nyguist 
Stability peint and indicating a system instability. Using 
equation 6.1the pole placement and robustness design method 
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Re 


Figure 6.4 Nyquist Diagram of 1/s+1. 


can be demonstrated . The first step is to establish a ktase- 
line for the design. The state unity feedback model in this 
Simple proklem gives pole locations of -2 and -2. This set 
of eigenvalues -2, -~2 were chosen as the baseline for the 
System. Since a seccnd order system reyuires only two feed- 
rkack gains to place the poles, the diagonal feedback gains 
were chosen for pole placement purposes. The pole placement 
and robustness design program was then used to place the 
Foles and in the process obtained the gains required to do 
this (1,1). A plot cf some of the Singular value criteria 
cbtained is shown in figure 6.6. Pole-zero plots of the 
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Figure 6.5 Perturbed Systen. 


closed- loop transfer functions of the closed-loop transfer 
matrix are shown in figure 6.7. In this case the feles and 
a zero are clustered about the -2 point and in the input two 
to output one channel a zero is located at the -1 foint on 
the pole-zero diagram. This point closely corresponds to the 
Rinimum Singular value frequency. The only significant 
observations are the ability of the pole placement and 
robustness routine to place the poles andthe relatively 
poor singular values indicative of low robustness. 

Since there is a reguirenment for two feedback parameters 
to set the poles of a second order system there should te no 
additional freedom in design to account for robustness. 
Case two was a run tc demonstrate this fact. Again allowing 
only two design variables for the problem the pole placement 
and robustness deSigr program waS run but with an objective 
functicn formulated to adjust the singular value above a 
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Figure 6.7 Closed-loop Poles for Pole Placement Only Case. 
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base level of 0.6. This Level would ccrrespond, using the 
universal gain margin chart, to a gain margin of -4 db to 8 
db and a phase margin of + 35 degrees. With only the two 
design variables to work with the pole placement and roktust- 


ness design program was unable to piace the poles and adjust 


the singular value level to the required value. The roles 
were placed at -2, -2 but the singular value minimum was 
Serr cn the order of -23 db or 0.07. The pole-zere plot 
remained aimost unchanged. ie Tse clear that additional 


degrees of freedom fcr the pole placement and robustness 
design rogram must be opened if robustness is to be 
accounted fcr. 

Case three was then run on the pole placerent and 
robustness program fry adding on additional degree of 
freedon. This case used the f£ gain as the additional 
design variable. A good choice as will be seen. Allowing 
the oftimizer routine the extra freedom to adjust the addi- 
tional feedrack gain term an excellent design was fcund. 
The singular value minimum LlLecame 0.885 as shown in figure 
6.8. Using the universal gain margin chart this corresponds 
to gain and phase margins of -6 db - 18 db and +52 degrees 
respectively. This is a considerabie improvement cver the 
original design. Tne factor that changed the design was the 
upper diagonal feedrack term which provides a cancelling 
factor for the cross-coupling term D yo = 50. This can be 


seen by looking at the system matrix eguation 6.11. 
= A- BE (6.11) 
where eguation 6.1z2 gives A. 


£9 ~I-2, -£,,-5,,f22 (6. 12) 
0 minis. 
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If the upper diagonal term of the matrix can be driven to 
zero the system will be decoupled into a diagonal systen. 
To place the poles cf this system at ~-2,-2 the diagonal 
feedback gains must [re one. The optimizer feedback gains 
for this case were f,, =1.00006, f,,=0.99998 and, most impor- 
tantly, f,, =-51.83401. Using these values in the A matrix 
A = -2.00006 1.83501 
0.0 - 1.99998 
The value in the upper right position in the system matrix 
has keen lowered considerably from the value of near 50 that 
appears in this position in the low singular value cases. 
lowering this system gain value decreases the cross-coupling 
perturbation effects cn the system. This can be more graphi- 
cally demonstrated by figures 6.9 and 6.10. PR tigume 6.9 
the transfer function shows a high gain of approximately 35 
db and a bandwidth of 50 rad/sec. In figure 6.10 this gain 
has been reduced to 6.0 db with a bandwidth of almost 2 
rad/sec. Two things are indicated by the figures, one, the 
open-loor Bcde plot of the transfer function of the cross- 
coupled channel can be used to indicate the robustness 
problem as evidenced Ey the high gain and bandwiith relative 
to the other feedback gains and, two, the mechanism used by 
the pole flacement and robustness design procedure to 
recover robustness is to reduce the relative gain and asso-~ 
ciated randwidth within the affected channel. Figures 6.11 
and 6.12 which are for the input one to output one channel 
are approximately the same as are the Bode plots’ for the 
input two to output two channel which are not shown. These 
figures indicate that no problem exists in the diagonal or 
direct ccourling terms. The pole~zero diagram for optim za- 
fon run case three is included in figure 6.13. The only 
Significant change in this flot as compared to figure 6.7 is 
the nevement of the zero in the off diagonal pole-zero flot. 
The zero is seen tc shift to one of the pole lccaticns 
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Figure 6.13 Closed-loop Pole-Zero Plot Case 3. 
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G@ueang the optimizaticn. PivcezZzeeOusiilrt bas tne effect of 
smoothing the frequency response curve in the vicinity of 
the frequency of the ainimunp singular value providing a more 
unifcrm gain distribution around this point. Table 1 shows 
comparison results fcr the feedback gains for s2veral cases 
of this kEasic problen. 

The logical extension of the problem to the general case 
is to allow all four feedback gains to become design vari- 


TABLE 1 
Comparative Results Simple Problem 


Feedback Gains 


variable £51> -51.83401 


variable £1 2=0 0.0 -0.01071 


: ee _ 
£54 variable i> SO. 50.6 0.00025 


All £'s variable - -60.4718 0.00280 | 





ables. The pole placement and robustness routine can then 
use full freedom in choosing all four of these feedback 
gains to ccmpensate for any cross-coupling effects within 
the system. Based on the previous analysis the two diagonal 
gains would be expected to approach 1 while the upper off 
diagonal gain moves tc -50 and the lower one moves to 0. The 
flexibility of the ALS program was reguired for the four 
design variable study. Several runs were made with various 
combinations of starting parameters for the feedback gains 


and oftizizer routines before a good design for the case 
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incorrorating all fcur feedback gains as variable was 
obtained. While this points out one of the linitaticns of 
optimization design routines the program was able to develop 
an improved design cver the baseline while employing four 
design variables. In the design produced for this case the 
Singular value level was placed above 0.6 for a gain margin 
of -4 db to 9 db and phase margin of + 35 degrees. In this 
formulation the optimizer was able to place the poles and 
meet the design Singular value level. The feedback gains 


produced by the optinizer were: 


f-1 1926 -60.4718 
0.0028 0.75478 
which are approaching the analytic design gains. After 


obtaining the feedback gains the OPTSYS program was used to 
obtain the necessary data todo a closed-loop pole-zere map. 
This plot is shown in figure 6.14 These plots indicate a 
Similar fpole-zero Ilccation to that found in the previous 
three design variable problem. Again the gain in the 
affected channel has reen reduced to compensate for the high 
cross-coupling perturkation within the systen. 

The design studies presented to this point have been 
based on breaking the system loop at the input as shewn in 
Ergure 6.15. In multivariable theory the location of the 
break in the loop changes the return difference fcr the 
system and the transfer function formulation. In the figure 
Number 1 depicts a system with an input loop break foint 
while number 2 is an output loop break point for output 
return difference determination. The return difference func- 
Beomefor the point 1 is written as I+FG while the return 
Beererence fcr point 2 is I+GF. The baseiine system not only 
has low singular values of the input return difference 
Matrix, the lowest being 0.0706, but also has low singular 
values for the output return difference matrix of point 2 of 


about the same order. To demonstrate the versatility of the 
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Figure 6.14 Pole-Zero Plot Simple Problem (4 Wariables). 
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Figure 6.15 System Block Diagran. 


pole placement and trcbustness routine a design was made in 
which tke output singular values were specified instead of 
the input values. Tke pole placement and robustness routine 
produced a design with gains that provided a highly robust 
system on the output side with gain and phase margin cf -6 
db to ‘and phase margin of +60 degrees. The design was not 
robust on the input side. Thus, designing for robustness at 
one point in the system does not necessarily give rotustness 
at all pcints within the systen. 

Cne final case that should be discussed is that of 
setting Foth input and output robustness criteria at the 
Same time. Excellent results were obtained for this case. 
The design routine flaced the poles at -1.97+.009j with 
feedback gains of 

0. 357393 -45.91757 

0.00425 Cecile? 
The input singular values were raised to a level of 0.74984 
or -4.5 db to 12 db gain margin and +43 degrees phase 
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Margin. The output values were raised to above 0.822 which 
corresfonds to -5 to 15 db gain margin and +49 degrees rhase 
margin. 

Tc Summarize, it can be stated that the rokustness 
problem for this system exists in the upper cross-courling 
channel (input two,output one). The lack of robustness can 
be discovered in two ways. The first method is to flet the 
Open-loor Bede plots cf each element of the transfer matrix 
and lock for extremely high gains and bandwidths relative to 
the cther transfer functions. The second method examines 
the singular values of the return difference matrix for 
Magnitude. Low Singular values correspond to low rokust- 
ness. The fole flacement and robustness deSign routine can 
increase rokustness Ey modifying feedback gains to reduce 
the effect of cross-coupling within the system. Observing 
the gain modification made by the pole placement and rotust- 
ness reutine the critical channel within the system that 
affects the robustness may be determined from the Bcde 
plots. The pole placement and robustness routine feedback 
gain changes also cause zero shifts during the rokustness 
Lecovery. The gain on the open loop Bode plot for the 
affected crcss-coupling channel is adjusted and the clcecséed- 
loop zeros, as seen cn the pole-zero diagram, are shifted. 
This zerce shift is ina direction which will combine with 
system poles to smooth the frequency response diagram in the 
Vicinity of the mininum singular value. 
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VII. A HELICOPTER STABILITY PROBLEN 


This chapter will’deal with a more practical application 
cf the numerical optimization progran. In this proklem the 
combined pole placement, robustness design proceiure will be 
applied to the linear lateral dynamic channels of a CH-47 
helicopter. The model is a highly coupled two-input two- 
cutput system that has been studied for its basics robustness 
Characteristics [Ref. 16}. The usual procedure for design 
of highly ccupled systems is to obtain a diagonally dominant 
closed-lcop system. This diagonal system will be stable but 
not rebust to cross-~feed parameters. Sandeil,et al, precduced 
three designs. Two cf these designs, while meeting Lasic 
performance criteria, had poor robustness. The third design 
was a relatively gocd design. The numerical optimization 
technique developed in the thesis was applied to the two 
poor designs and shown to provide substantial improvement in 
robustness. 

The systems were designed to satisfy specifications to 
step input response and stability margins as stated in mili- 
tary specifications. Specific design parameters for each of 
the three designs presented were not available. The designs 
were all ccnsidered to meet the performance specification 
criteria and stability margin freguirements. It was shown 
that two of these designs were extremely sensitive te modei 
errors. The classic Nyquist technigues did not predict this 
sensitivity. The singular value analysis did indicate sensi- 
tivity problems. 

The model is that of a CH-47B helicopter lateral dynamic 


system in hover. The dynamic nodel of the system is 
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ine 
T] 
> tps 
b 
+ 
eu 
SC 


(7.1) 


X =(V,p,I,¥) (jez) 
a = (Sg0 5c) Cleo) 
where 
Lae | -1.42 19/54 5) 39.2 
A = 0.01 -0.7 -0.07 0 
0.04 ~0.05 a0 275 0 
| 0 1 O. 11 0 J 
and where 
Oe 12 0.95 
B= weou —8.37 
Ae 
fou 0.02 
0 0 


with full state available for feedback. Table 3 is a summary 
of parameters. The system is not open-loop stable. 

Three ccntrol laws are formulated to satisfy the desired 
perfcrmance specifications. Equation 7.4 is the Easic 
contrcl law. 


= ~F x i hy Pe (7.4) 
m=1,2 0r 3 

where the following values of F and h were are: 

F = He 5 70..6 —- 


i 
0.024 -2.71 02266 0-7 ae 


E> = | 0.198 184.0 18.3 142.0 
-0.01 -1.592 -0.189 -1.47 
F, = 0 0 25.5 0 
0 -4 0 -27 





TABLE 2 
CH-46 Helicopter Parameter Definitions 


ft/sec} Vehicle body y-axis earth relative velocity component 
Pe G/sec Roll rate 
rad/sec | Yaw rate 

rad Roll attitude angle 

in Yaw rate rotor deflection control 


in Roll rate rotor deflection control 
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The Ye is a step input command that must track. Figure 
7.1 is a diagram of the control structure for the feedback 
contrel laws. Figure 7.2 is a detailed layout of the 
systen. All three cf the designs have negative real eigen- 
values and frovide stable overdamped responses. All three of 
the above designs meet the desired design specifications. 
The loop Bode flots, [Ref. 16], indicate all three 
designs to have good stability margins when considered a 


loop at a time. Since these designs were all full state LQ 
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Figure 7.1 Feedback Control Structure. 


as 








System Diagragm. 


Figure 7.2 
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designs with a diagonal control weighting matrix they should 
posses at least -6 dr to infinite gain margin ani 60 degrees 
of hase margin. When the singular values of the taree 
designs are computed a robustness problem is indicated by 
low singular values cf two of the designs. Figure 7.3 pres- 
ents singular value flots of all three designs. This plot 
Shows that designs 1 and 2 both have very low minimnun 
Singular values for the return difference matrix. Design 1 
goes as low as -20 db near 10 rad/fsec in freyuency while 
design 2is down to -34 db at ‘freguencies up to 1000 
rad/sec. Design 3 is a good design with singular values that 
remain above one thrcughout the freguency range of interest. 
Using the universal gain and phase diagram as discussed in 
Chapter 6 this equates to a gain margin of about -6 db to 
infinity anda phase gzargin of 60 degrees as expected. 

The system stability in design 1 may be affected by 
perturfaticns occurring in the actuators as shown in figure 
7.4. iieene Output axis coupling from dc Sprlis into Jp in 
the frequency range from 0.5 to 50.9 rad/sec with a magni- 
tude cf 0.12 and a fhase of 60 degrees then the system can 
become unstable. This could be caused by nonlinear terns, 
been parts, or system saturaticn. | 

The second design as shown in figure 7.5 may have a 
stability derivative variation between Sg and p. If this 
derivative varies frem about 0.04 to -0.96 the system can 
become unstable. The center of gravity location, trim of the 
gic ract and rotecr coupling can all affect this lerivative. 

The pole placement and robustness design technigue 
developed for this thesis research was applied to designs 1 
and 2 tc oktain an improved robustness for these designs. 
Design 1 will be considered first. In this design a 
cross-feed perturbation through the actuator can prcecduce 
instability as indicated by the low singular values at lower 
frequencies. To study this problem with the pole placement 
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Figure 7.5 
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and robustness design method it is assumed that the eigenva- 


lues (pele locations) of the system as developed in 
fRef. 16}. are the reguired poles for the performance 
@m@a teria. Once the pole locations are set, rotustneéss 


criteria must be selected. From the universal gain and fhase 
Margin curve discussed earlier two choices of singular value 
levels were made for this problem. The first singular value 
level chosen was 0.6. This corresponds toa gain margin of 
-4.0 db to 8 db and a phase margin of about 35 degrees. The 
second value chosen was 1.0. With corresponding gain margin 
of -€ db toc® and a phase margin of 60 degrees which are 
the characteristics of aLQ regulator design with diaconal 
weighting matrices. 

For singular value level 0.6 the pole placement and 
robustness design routine places the poles as Shown in table 
2 The slight differences in these pole locations appear to 
have insignificant effect on the performance as Snown ir the 
response curves of the system. The feedback gain adjustment 
moves the minimum Singular value from about 0.11 with very 
poor phase and gain zargins toa level of 0.66. This is 
above the desired vaiues of gain and phase. The imprevement 
in rekustness came from modifying the feedback gains in 
channel pr. By greatly reducing the gains in channel dg the 
optimizer minimizes the influence of the cross-coupling from 
channel $c. In this way a much larger spill over of channel 
bc may be tolerated through the actuator before the systen 
will keccue unstable. 

The primary mechanism of robustness improvement in this 
Peopsecm was a reduction in the gain levels of the affected 
channel. The feedback gains presented in table 4 show the 
modificaticn of these gains from those Hite Zed. 1h 
[Ref. 16}. While all the gains are modified the f,g gain 
undergoes a significantly larger change than the cther 
design 1 case 1 feedkack gains. Looking at the open-loop 





TABLE 3 
Design One Pole Placement 


Desired Location Actual Location 


1 mate say 7 BO Tisi 8) S 












mit. 3O30 -12.0083 
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transfer furctions cf this optimized problem and comparing 
them with a non-optimized set of transfer functions for 08 
to as shown in figure 7.6 and 7.7 show very iittle 
change in the system gain, both are around 42.0 db, anda 
bandwidth of abcut 20 rad/sec. There is a significant 
change in the phase diagram which is caused by the zero 
Meeation shift. 

The transfer function for Oc to de depicted a gain 
increase of about 3 db for the optimized design while the 
phase remained similar for both transfer functions. In 
transfer function dc to dg some important aspects of the 
problem are observed. The Bode diagram of the open-loop 
transfer function of de to Sg clearly indicates the 
cross-coupling problems and the pole placement and rokust- 
ness design routine's mechanism of optimizing the systen 
gains to increase robustness. The gain is reduced frem 93 
db to akout 78 db and the bandwidth is reduced from above 
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100 rad/sec to about 30 rad/sec. The reduction of bandwidth 
and gain in the loop be to OB yields an increased toler- 
ance to perturbation. The transfer function for Ze to 
Og is shown in fires! /.oeaud 7.9. in the transfer “iunc- 
tions for 582 to Se the bandwidth is slightly increased 
from 0.6 to 0.8 rad/sec and the gain actually increased fron 
cede to 14 db. The big change in the overall systen, 
however, is in the transfer function from Se to oa . This 
is the channel that the destabilizing perturbation enters 
and by greatly reducing the gain and bandwidth in this 
channel through a change in feedback gain, the optimizer 
routine has brought the entire system gains to more Lalanced - 
conditiors and recovered a highly robust design. 

The gain changes associated with the robustness improve- 
ment cause the zeros of the various closed-loop pole-zero 
diagram tc move. A comparison of the eight pcle-zero 
Giagrams is shown in figures 7.10 to 7.13. The significant 
feature of these pcle-zero diagrams is the shift of the 
zeros of the optimized design ina direction that attempts 
to equalize or balance the frequency response for freguen- 
cies in the vicinity of the Pinimum singular values. the 
pole-zero diagram of SB to v will be discussed as an example 
of this effect. In 7.10 the nonoptimized zeros are located 
Beeut —2 + 4.5j and -19.8. When the pole placement = and 
robustness routine has completed the feedback gain mcdifica- 
tion these zeros have shifted near <-11and -6+3j. =the 
effect of these zero shifts is to combine with the pcle 
locations tc equalize the freguency response as depicted in 
figures 7.14 and 7.1&. Zero shifts for the remainder of the 
transfer functions provide similar results in the cther 
channels. By moving toward the frequencies associated with 
the minimum Singular values the zeros have balanced the 
overall frequency response of the system in each channel. 


While the channel gain modification is the primary mechanisz 
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Figure 7.13 Pole-zero Plots for Design 1 (cont.). 
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Bere LObustness recovery the zero shift associated with the 
feedkack gain changes is directly related to the overall 
freyguency response of the systen. 

In design i case 2 the minimum Singular value level was 
increased from 0.6 to 1.0. The pole placement and rokust- 


ness design procedure also gave a better design for case 2 


TABLE 4 
Helicopter Problesa Feedback Gains 









f 28 1p at 
GAIN VALUES Be gle fl elt 


fees 
APR ee 24 
One Stemi y7ee 2615656 77 96629 ~32.91595 
Case 1 -0.00567 -2.55646 0.39039 -15.04805 


-6.36249 -1.53746 74.20387 
0.00813 -2.60073 0.65395 -13.47128 


















-0.96660 -4.54597 0.75436 11.96808 
0.00916 -3.03069 -0.36391 -0.34991 








Ey modifying the feedback gains as shown in table 4. The 


Change in the transfer functions for Oe to OB (figure 
7.16 and 7.17) is similar to the change seen in figures/?7.8 
ona 7.9. The gain was reduced from above 90 db to about 75 
db and the ktandwidth cut from above 100 rad/sec to about 25 
rad/sec. The added requirement of increased robustness did 
not result in a Significant change in performance. 
{Ref. 16]. 
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Tesign two which is the design that is subjest to paran- 
eter uncertainty in the dg to p channel was also studied 
using the pole placement and robustness design routine. In 
this design the singular values were very low until ireguen- 
cies akove 100 rad/sec were reached. To improve the design 
it was assumed that the pole locations that corresponded to 
{[Ref. 16] were the reguired pole locations and that a 
robustness singular value level of 0.6 would provide 
adequate gain and phase margin for the design. The pole 
placement and robustness routine adjusted the gains in this 
problem until a minimum singular value of 0.6 was oktained. 
During this adjustment the gains in channel Jg were consid- 
erably reduced to offset the cross-coupling between the two 
channels. Again, in this design as before the channel 
cross-coupling between de and OF shows a marked change 
in bandwidth and gain from above 100 to 35 rad/sec and about 
78 db to 22 db respectively. Figures 7.18 and 7.19 shew the 


Meansicer function plets for this tera. Figure 7.20 shows 
the singular value improvement. In figure 7.21 the time 
response is plotted. As can be seen in the plot the 


improvement in robustness for this problem results in very 
Sluggish response and degraded performance. 

To summarize for this problem, a given performance level 
has been chosen in terms of pole locations. The level of 
robustness has been set for a desired gain and phase margin 
kased on the universal gain and phase margin curve. The pcle 
placement and robustness routine has been able to improve 
the robustness level by changing the feedback gains that 
affect the chamel fe GEGss- Coupling. This rokustness 
recovery is affected by modification of the system feedback 
gains in such a manner that cross coupling gains are reduced 
so that small cross-ccupling perturbations do not drive the 
system inte instability. The open-loop transfer function 
plots have keen used to indicate how this mechanism operates 
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Figure 7.21 Time Response Design Two. 


111 


40.0 





and have keen shown to be an alternative indicator of chan- 
nels that may be affected by cross-feed perturbations. The 
Pole-zere diagrams. of the closed-loop transfer functicns of 
the transfer matrix further indicate that zero movement is 
in a direction that egualizes the gain level of the 
frequency response curves in the vicinity of the lowest 


Singular values providing a more balanced system response. 





VIII. SIMPLE OBSERVER 


ae = ab ee 


The ¢ole placement and robustness design procedure can 


also be used for robustness recovery in observer design. 


Say Nersaa 


aC xX 





A 


X=AX+KY + BU 


Figure 8.1 Simple Observer. 


Given an observer as represented in figure 8.1 it has been 
shown that the system differential eguation may he written 


as eguation 8.1 


x =/[A -BFi/fx) + fa} cr (8.1) 

3 .. A.-B i R ;| 
A. = A - KC 

where x is the state and x is the observer variable, and the 

solution may be separated to independent solutions for the 


feedkack gains and the observer gains. This separation 
allows the use of the feedback gains to set the pole 
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locations and then a second optimization run with set feed- 
rack gains so that the observer gains may be computed to 
adjust tke system rokustness level.The observer pole loca- 
tion could also be flaced using the observer gains, kK. In 
the currently implementation of the pole placement and 
robustness design reutine the observer poies are simply 
restricted to areas cf the left half plane. 


In this chapter a simple stable observer system will be 


analyzed based oon [ Ref. 17]. Given the system, equation 
8.2, 
_¢ =| (OSS > ai 0 So ay (8.2) 
=5) ah 1 -61 


y= [2 1] ey 
where EF (¥)=E(4)=0 and E(¥(2)y (et) ) =E (4 (4) 7 (e)) = FS (t-7) with 
the feedtack law of equation 8.3 


u = -(5C 10)x + 80r (8.3) 


An analysis has been done to ccmpare results of the numer- 
ical optimization procedure with results presented in 
fRef. 17}. The results for an optimal regulator design 
uSing guadratic cost criteria as stated in equation 8.4 are 


given in table 5 
J = ftatax uz) dt (8.4) 


with p= 45 (735 1) 

Figure &.2 shows a Nyquist plot of the full stat2 regulator, 
the optimal filter and a fast filter. The full state design 
had poles at s= -7.0+j 2.0 and feedback gains of 50 and 10. 
mre, OEtinal filter as shown in table 5 had poles of -7.0 + j 
2-0 with gain and phase margins of ~6.75 db and +15 degrees. 
irensortimum filter gains were 30 and -50. Using a faster 
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Figure 8.2 Nyguist Plot. 


filter also gives pcecr gain and phase margins. The gain 
Margins are on the order of -.98 while the phase margin is 
less than 10 degrees. The bandwidth also increased fren 12 
to 40 rad/sec. A reccvery procedure based on a modification 
to the process noise matrix [Ref. 17] may be applied to the 
Froblen. This procedure can recover a large amount of the 
robustness that was lIcst with the observer addition. Figure 
&€.3 shows data obtained for several trials of thea fictitious 
noise procedure. The gain, phase margin and cthker 
parameters may be found in table 5. 
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TABLE 5 
Observer Parameter Data 


FILTER GAIN PHASE ERROR STATE FILTER 
POLES MARGIN MARGIN COVARIANCE COVARIANCE BAjh 
db deq E [tx-s)(u-n)"] wen 


Optima} Q7 “163 221 -613 30 
LQG Design 27423 15 
-163 277 | -613 2070 . 


Fast Filter Adjust. 
ment Procedure 
















6280 
-12250 


-12200 
23800 















Fictitious Moise 
Adjustnent Procedure 


q = 100 





°570 -613 
-570 1170 -613 


POPLAR CASE I 1.00018 GM=-5.5 dB 
Elie PM=117 deg 
POPLAR CASE II 18,22148 GM=-19 dB 
32.59476 PM=58 deg 


The singular values of the state feedback and optizal 
cbserver systems were computed for comparison s#ith results 


produced using the role placement and robustness design 


recovery procedure. Figure 8.4 indicates a loss in rckust- 
hess represented by the observer singular values. Lower 
Singular values are less robust. Por this single-input 


Single-outfut observer the Nyguist diagram will be used to 
define gain and phase margins. Since a singular value of 1.0 
is indicative of a linear quadratic level of gain and fhase 
Margin,i.e. GM=-6 db, and PM=460 degrees, this was chosen 
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Figure 8.3 Hyguist Plot for Robustness Rerovery. 


as the design level. The pole placement and rokustness 
routine was used to recover robustness while setting the 
pole placement at the state feedback pole locations of 
-7ij2. The minimum input singular value level was set at 
1.90. The pole placement and robustness routine was also set 
to place the observer poles anywhere between -100 and -2 
that would provide robust design. The plant poles were 
placed at -7.05 +j 1.82. The feedback gains for this run 
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were £0 and 10.09. The optimizer output prodused a filter 
With a gains of 1.00018 and ~2.02742. The observer pcle 
locations were -1.99 + j .001, very near the plant zero 
location. A Significant change was produced in singular 
values. Figure 8.5 shows all singular value curves plctted 
together. The optimizer solution for this problem is well 
above the optimal filter curve at low frequency. Figure 8.6 
shows the Nyguist plct of the optimizer developed design. 
The system has a gain margin of -6 db and a phase margin of 
117 degrees. The mest significant differences between the 
two desicns being that the observer poles are close to the 
plant zero locations and the filter gains are much lower for 
the cptimizer soluticn. 

Using the OPTSYS program with the pole placement and 
robustness routine ccmputed gains as the design parameters 
the data for the observer filter was computed. The error 
covariance matrix was found to be 

BOG o ~456.5 

-456.5 700.1 
These values compare favorably with the trends established 
Emeeacle 5. The last comparison of the pole placement and 
robustness design reccvery procedure was an analysis of the 
time retrcnse curves. Figure 8.7 shows the comparison flot. 
The design obtained using the pole placement and roktustness 
routine did not degrade system performance. 

Cneé additional analysis was conducted that set the 
desired design parameters slightly differently. In this run 
the peje placement and robustness routine was set to plece 
the okserver poles between -10 and -100. The pole placement 
and robustness routine was unable to totally satisfy this 
requirement. It violated one of the constraints and noved 
the larger cbserver role to -2.25 which is near the orftimun 
pole locaticn and also the high gq? values of the fictitious 
noise precedure. The smaller pole was moved out to -70.8 


119 





0 


Ww 
FULL STATE 
© 
a POPLAR 
OMS 10h 

a \ 

>) a 
7 
=] ~~ 
«wt! “XN 
> \ \ 
{Oo POPLAR ‘ 
Sa CASE I ‘ 
S \ 
ips \ 

i 
\ \ 
‘\ 
© aa ae — 44g Eee 
q OPTIMAL . - Sa 
OBSERVER 
© 
0 LAG 100.0 


10.0 
FREQUENCY = RAD / SEC 
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Figure 8.6 Nyguist for Cosputed Robustness Recovery. 


which corresponds to the same level of pole movement found 
when g@* reccmes large. The singular value levels were raised 
Significantly as shown in figure 8.5. The state and error 
covariance matrices were of the same order as the matrices 
found in table 5. This design has a phase margin cf 58 
degrees. The gain margin was about -19 db. Even though scne 
of the ccnstraint conditions on the design were not met the 
design demonstrates excellent robustness recovery. 
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Figure 8.7 Time Response Plot for Simple Observer. 
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The pcie placement and robustness design routine has 
produced a robust design ‘for this observer based systen. 
This design is oktained by uSing a numerical optimizaticn to 
directly Tanipulate the feedback and filter gains. 
Modificaticn of the 1Q functional equation as done in the 
fictitious noise procedure is not reyuired. The pole rlace- 
Ment and robustness routine solution for this problem has 
lower gains than those found using the fictitious noise 
adjustment. The routine prevides a good, direct methodology 
for selecting the feedback and filter gains for a rcbust 


observer design with excellent performance. 
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This chapter will be devoted toa short discussicn of 
the rerustness recovery of a fourth order observer tased 
contrcller. The problen is a helicopter problem, [{Ref. 6]. 
In this case the helicopter model is that of the longitu- 
dinal ccntrel loop of a CH-47. The nominal model is taken to 
be the system of equations 9.1 and 9.2 for an aircraft speed 


Smeecrty Kncts. 


fe} -0.02 0.005 2.4 -32. |x +}0.14 -.12 ]u (9.1) 
"0.14 0.44 -1.3 -30. 36 -8.6 
Os 0.018 -1.6 1.2 5m 002 
On 0. ike 0. Ca Oe 
J 
me= 10 10 0 x (S752) 
io 0 0 * 


in this prceblem the controller is formulated as in figure 
9.1 which leads to an open-loop transfer function of the 
form of equation 9.3 


K (S$) G(s) =F (SI-A+ BF+KC) 71 KC (SI-A) -1B (9.3) 


The pcle placement and robustness recovery procedure was 
eppiiéed to this fproblen. 

Pirst, the standard full state feedback design was 
carried out using the Naval Pecstgraduate version of OPTSYS. 
This design produced excellent singular value output for the 
return difference as shown in figure 9.2. The lewest 
Singular value being essentially 1, corresponding to a LQ 
design with -6 db to @ gain margin and 60 degrees of fhase 
margin. The time resrcnse of the system was good as shown in 


ee 








Figure 9.1 Observer Based Controller. 


figure 9.3, reaching steady-state in about four seconds with 
only a slight overshcect. Since the characteristics for this 
system are acceptable no further design iteraticn was 
carried out. The full state feedback became the baseline 
design. Assuming that full state feedback was not availalle 
and only two measurements could be produced an observer was 
developed to control the two measured outputs, vertical 
velocity and pitch attitude. USing the neasurement matrix of 
equation 9.2, OPISYS was used to develop an optimal cbserver 
system for this yproblen. The Singular values of the return 
difference matrix and the time response were plotted for 
compariscn to the full state design. These are shown in 
figures 9.2 and 9.3 respectively. Note that the singular 
value of the return difference matrix is as low as 0.16 at 4 
rad/sec. This equates to a gain margin of -1 db to 1 dk and 
less than 10 degrees phase margin. The time response is 
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Figure 9.2 Singular Value Plot of Observer Results. 
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Flotted in figure 9.3. The steady-state was reached for the 
velocity after a slight cvershoot in about four seconds. 

As the final step in the analysis the pole placement and 
robustness routine was employed to recover robustness cf the 
observer based system. The pole placement and robustness 
routine was first used to set the poles at approximately the 
same location as the foles of the LQ regulator. These were 
assumed to re the desired pole locations. The pole placement 
and robustness routine was then used to vary the filter 
gains, K, until the desired level of robustness was reached. 
The desired singular value level chosen was 1.0 which ccrre- 
sponds tc -6 db toec®O db gain margin and 60 degrees of fphase 
Hargin. Figure 9.2 shows that the pole placement and rotrust- 
ness procedure failed to totally recover the rokustness 
level to 1.0. The minimum singular value reacaed was only 
0.936. This singular value eguates to a gain margin of -5.0 
db to om anda phase margin of 55 degrees. While this is 
slightly less than the design objective, it is far superior 
to the optimal observer design discussed inthe frevious 
paragraphs. This recovery was made by making an optiniza- 
Piomerun, finding the filter gains with near zero values and 
freezing these values at zero. This reduced the number of 
design variables the pole placement and robustness routine 
waS reguired to manifulate in a second optimization run and 
gave a higher robustness solution. 

This analysis and the second order observer analysis 
presented earlier clearly indicate that it is possible to 
use the pcole placement and robustness procedure and separa- 
tion principle to develop a robustness recovery procedure. 
The pole placement and robustness routine provides rokust- 
ness reccvery by direct modification of the faedback and 
filter gains. This procedure requires no nodificaticn to LQ 
cost functicnals or cther parameters as done in the ficti- 


tious noise adjustment method commonly used for rotustness 
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CTeCOVELy. By providing direct gain adjustment the pcle 
placement and robustness procedure results in a practical 
design with relatively low observer gains and good ;ferforn- 
ance. The procedure is Simple and straight forward with the 
only difficulty being a reguirement to sometimes nmcdify 
initial starting values or optimizer codes to force the 


soluticn toward the desired point. 


tS, 


X. CONCLUSIONS 

An effective method of roktustness multivariable control 
design utilizing a numerical optimization based algorithn 
has Eeen developed. The pole placement and robustness design 
routine courfled with the Automated Design Synthesis program 
provides the designer an excellent tool with which te attack 
the rebust design prctlen. 

The pele placement and robustness deSign coutine has 
demonstrated the capability of providing designs that solve 
the problems caused by cross-coupling perturbations which 
reduce rokustness in multivariable systems. This design 
improvement is accomplisned by modifying the system feedback 
gains in such a manner that the gain in channels that are 
affected by cross-coupling perturbations is egualized witn 
other system gains to reduce this cross-coupling effect. The 
gain changes are accompanied by zero shifts which also 
influence the gain distribution and freguency respense of 
the systen. 

Ferturbation proklems in multivariable systems have been 
Shown to be detectakle by singular value analysis and by 
using the Bode magnitude diagram of the open-loop transier 
functicns of the system. In the open loop transfer function 
large differentials in Bode gains and bandwidths are indica- 
tive of problem areas for cross-coupling perturkations. 
Robustness is obtained by the pole placement and roctustness 
design program by modifying those gains and bandwidths asso- 
ciated with the cross-coupling perturbations thus reducing 
the amount of energy coupled from the perturbation into 
other channels. An associated zero shiit has been observed 
when these gain modifications take place. This zero shift is 


in the direction of roles that are located in the vicinity 
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cf the freguency of the minimum singular value and tends to 
equalize the frequency response curve gains in this region. 
The use of numerical optimization to recover rotustness 
in okserver based designs was demonstrated. The pole flace-. 
ment and robustness routine was applied woe CEOLLlens 
previously solved using the fictitious noise procedure for 
robustness recovery. The direct manipulation of feedback and 
filter gains by the fole placement and robustness routine 
rrovided a highly rokust design with relatively lew filter 
gains. The problem of robustness BecOVery @rin fi11ter- 
ebserver designs has been soived in a straight forward and 


highly practical manner. 
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